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We develop a framework to study the phase space of a system consisting of a scalar field rolling 
down an arbitrary potential with varying slope and a background fluid, in a cosmological setting. We 
give analytical approximate solutions of the field evolution and discuss applications of its features 
to the issues of quintessence, moduli stabilisation and quintessential inflation. 
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I. INTRODUCTION 

Scalar fields play a central role both in particle physics 
and cosmology. For example, in the Standard Model 
of particle physics, the Higgs boson generates particle 
masses through a mechanism of symmetry breaking. Su- 
persymmetric models, believed to solve the hierarchy 
problem and to suggest a grand unification scale, add 
a whole new set of scalar (and fermionic) partners to 
the particles of the Standard Model. In string theories, 
the only physical constant is the string tension being all 
the other constants generated dynamically through scalar 
fields, the moduli. 

In cosmology, a scalar field, the inflaton, has been sug- 
gested to be responsible for an early inflationary period 
in the history of the universe. This scenario can account 
for the extreme flatness and homogeneity of our Uni- 
verse W . More recently, measurements of the apparent 
magnitude-redshift relation of type la supernovae (Snla) 
support indications given by the combination of cosmic 
microwave background radiation (CMB), galaxy clusters 
and light elements abundances measurements, that the 
Universe is presently undergoing a period of accelerated 
expansion H. Once again, the existence of a scalar field 
rolling down a potential has been suggested as an expla- 
nation for this late inflationary period. The scalar field 
is usually called "quintessence" pj. However, whether 
the inflaton and quintessence could be the same field, 
remains an interesting open question. 

Considering these and other examples of applications 
of dynamical scalar fields in physics, it seems important 
to take a closer look to the features of their potentials. 

For instance, in the issue of quintessence, if the poten- 
tials have attractor solutions {i.e. the late time dynamics 
is independent of the initial conditions) then there is a 
chance to weaken the fine tuning problem associated with 
a cosmological constant term in Einstein's equations pi . 
In the same way, when allowing for dynamical evolution 
of the moduli fields, the existence of attractor solutions 
opens up a larger region of initial conditions for which 
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the fields can have successful stabilisation at their vac- 
uum expectation value ||,|| ■ 

The aim of this work is to extend the results of |7|-|l0[| to 
more generic type of potentials, motivated by the issues 
of quintessence, moduli stabilisation and quintessential 
inflation and to give analytical support to some of the 
conclusions in [ pl| . 



II. SETUP 

We consider a spatially-flat Friedmann-Robertson- 
Walker Universe containing a scalar field <fi with poten- 
tial V(<j>), and a barotropic fluid with equation of state 
Pb = (7 — l)ps, where 7 is a constant (e.g. 7 = 4/3 for 
radiation and 7 = 1 for matter). The governing equa- 
tions of motion are, 



H = 
Pb = -3yHp B 
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subject to the Friedmann constraint 
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where n 2 = 8irG and dots denote derivatives with re- 
spect to time. The energy density and pressure of a ho- 
mogeneous scalar field are given by p^, — <j) 2 /2 + V and 

V<i> = </ ,2 /2 — V, respectively. 

Following j|], we define the variables, 
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where a prime denotes a derivative with respect to the 
logarithm of the scale factor a, N = In a. 

The effective equation of state for the scalar field at 
any point yields, 



1<P 



_ P <p+P<p _ T _ 2z 
P0 ~~ 4> 2 /2 + V ~ x2 +v 2 



(0) 



constrained between < 7^ < 2. In terms of these new 
variables the equations of motion read: 



-3s + \J \y 2 + \x [2x 2 + 7 (1 -x 2 - y 2 )] , 
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where we have defined 
A 
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Sec jllj and respectively, where these definitions were 
first introduced. Note that both A and T are, in general, 
<fr (and thus, time) dependent. 

The contribution of the scalar field to the total energy 
density, Q<p = n 2 p^jZH 2 = x 2 + y 2 is bounded, < 
Q<P — 1, if P > 0. Hence, the evolution of the system 
is completely described by trajectories within the unit 
circle. Moreover, since the system is symmetric under the 
reflection (x,y) — > (x, —y) and time reversal t — > —t, we 
only consider the upper half disc, y > in what follows. 

For a generic scalar potential, one can identify up to 
five regions in the phase space diagram (x,y). As an ex- 
ample, in Fig. |lj and Fig. |2| we show an exact numerical 
solution of the evolution of the system for a double ex- 
ponential potential. In these figures, region 1 represents 
a regime in which the potential energy rapidly converts 
into kinetic energy; in region 2, the kinetic energy is the 
dominant contribution to the total energy density of the 
scalar field ( "kination" ) ; in region 3, the field remains 
nearly constant until the attractor solution is reached 
("frozen field"); in region 4 the field evolves into the at- 
tractor solution, where the ratio of the kinetic to poten- 
tial energy is a constant or slowly varying; and in region 5 
the potential energy becomes important, the scalar field 
dominates and drives the dynamics of the Universe. 

In section III we will briefly discuss regions 1, 2, and 
3. However, in this paper we will be mostly concerned 
with regions 4 and 5 of the evolution since (as it will be 
shown) they correspond to stable solutions yet not having 
true critical points associated with them. As we will see, 
this feature plays an important role in many cosmological 
phenomena. 




FIG. 1. Evolution of the scalar field energy density p^, in 
a radiation background fluid ps for a double exponential po- 
tential. Regions 2 to 5 represent respectively, kination, frozen 
field, evolution in the attractor and scalar field domination. 




FIG. 2. The scalar field dynamics in the phase space di- 
agram for a double exponential potential. Regions 1 to 5 are 
the same as in Fig. h. The dashed line represents the unit 
circle x 2 + y 2 — 1. 



the A = 00 surface can be deduced by continuity from the 
A = 00 solution, even though the latter is not physical. 

We bring the plane A = 00 to a finite distance from 
the A = plane by the transformation 



A 
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0<5< 1 , 
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so that the 6 = surface corresponds to the A = surface 
and the 6 = 1 surface corresponds to A = 00. On the 
plane 6 = 1, we find 



III. EARLY EVOLUTION 

For a wide range of potentials (including the potential 
in Figs, [j] and H) , when cf> is high up the potential slope 
at the beginning of the cosmological evolution, then A is 
large initially. The behaviour of the trajectory close to 
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where r is a new time coordinate defined by dr = XdN 
Solving Eqs. QlCj), we obtain, 



= A tanh 



y = A sech 



- 2 (r-r ) 



(r-ro) 



(11) 
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where A and tq are arbitrary constants. 

As shown in Fig. 0, the trajectory for region 1 is nearly 
circular, well described by Eqs. ( |TT| ) and (|l2|). For A very 
large, the scalar field potential energy turns into kinetic 
energy very rapidly. 

In region 2, the trajectory is very close to the x-axis 
where the scalar field energy density is kinetic energy 
dominated (hence y w 0), and its behaviour in the x 
direction is described by 



i' = --(2-7)3(1 



x 2 ) 



(13) 



For x < 0, x' is always positive and x increases; for x > 0, 
x' is always negative and x decreases bringing the trajec- 
tory towards the A - axis (in the 3-dimensional phase 
space {x,y, A)). 

Close to the A-axis, the system is described by 



;(2-7)s. 
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Therefore, the trajectory in region 3 evolves along the x 



-3(2-7)^ Since the back _ 
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direction at a rate of x c 

ground fluid is dominant at this stage, H 2 ex a 
we have the kinetic energy falling off rapidly as 4> 2 ~ 
x 2 H 2 oc a~ 6 . The trajectory moves in the y direction as 
y oc a 37 / 2 , that is V ~ y 2 H 2 rs constant and the scalar 
field is essentially "frozen" . 



IV. THE TRACKER SOLUTION 

In this section we study the case in which the scalar 
field evolution is well approximated by a linear relation 
between x and y, as it happens along region 4 of the 
example illustrated in Fig. ^. In order to study the evo- 
lution of the scalar field when the slope of the potential 
is varying, i.e. A is changing with time, we rewrite Eqs. 
([?]) by defining e = 1/A, x = eX, and y = eY . 

We then have, 

X' = \/6(r - l)X 2 -3X + J |r 2 + 

-X[2e 2 X 2 + 1 (l-e 2 X 2 -e 2 Y 2 )}, 



(15) 



Y' = V6(r - 1)XY - J -XY + 

-Y[2e 2 X 2 + 7 (1 - e 2 X 2 - e 2 Y 2 )} , 
e = V6e(r - l)X . 



For e small (A large) or T rs 1, e becomes nearly constant. 
Moreover, if Y is also nearly constant we can solve X' — 
Y' = and find the "instantaneous critical points" , 



*o(A) = W^ 
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y 2 (A) = ^(2- 70 ), (16) 



where the equation of state of the scalar field is 



1, 



70 = -[ 7 +(2r-l)A 2 /3]± 



^[-7 + (2r - l)A 2 /3] 2 + 8 7 (r - 1)AV3 . (17) 

The plus root leads to unphysical results, so only 
the negative one will be used in what follows. The 
contribution of the field to the total energy density is 

^ - 3 70 /A 2 . 

Note that in the limit when r — 1 w one has 



70 = 7 



2(T-1) 



1 - 3 7 /A 2 



(18) 



In other words, for potentials with small curvature, the 
equation of state of the scalar field is very close to the 
equation of state of the background fluid, and it is said 
that the field "tracks" the background fluid. This expres- 
sion can account for values of f^ > 1/2. 

From Eq. (J17J), it can also be shown that in the limit 
of large A, when the background fluid is completely dom- 
inating, we recover the expression in B, 



70 



2r- 1 



(19) 



In the limit when T — 1 ~ 0, the above expression is 
equivalent to Eq. (^8|) in the limit of large A. 

One can now study the stability of the "instantaneous 
critical points" linearising Eqs. (|15| ) i n u an d v with 
X = X c + u and Y = Y c + v resulting in the first-order 
equations of motion, 



M 



(20) 



Using Eq. ( |l7J ) to write V — 1 in terms of 7^,, we find the 
eigenvalues, 



m± = --2 [(7 - 70)(370 + A 2 ) + (2 - 1<f> )\ 2 ] x 
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A 2 (2-70)(7A 2 ~3 7 |) 
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[(7-70)(370 + A 2 ) + (2- 7 0)A 



(21) 



It is straightforward to check that Eq. (|2l| ) reduces to 
Eq. (A8) of when 7 = 7^ (i.e. the pure exponential 
case) and to Eq. (15) of Q in the limit of large A. 

The system is stable if the real part of both eigenvalues 
is negative. From Eq. (El]), this is completely assured if, 



70 < 



7 A 2 



and 



, 3 7 ~2A 2 / 
T *< 6 I 1 



12A 2 (2+ 7 ) 
(3 7 -2A 2 ) 2 
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Moreover, if the quantity under the square root of Eq. 
( pl| ) is negative, the critical points are a stable spiral; 
and a stable node, otherwise. In Fig. || and Fig. |j one 
can see the dependence of the eigenvalues on 7 ^ and A 
for a background fluid of matter ( 7 = 1) and radiation 
( 7 = 4/3), respectively. The lines split when the quan- 
tity under the square root of Eq. ( |2l|) becomes positive. 
The first of the conditions in Eqs. (p2h marks the point 
for which only one of the real parts of the eigenvalues be- 
comes positive, and the second condition the point where 
both real parts of the eigenvalues become positive. In the 
limit of large A the latter is 7 = 1 + 7 /2. 
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FIG. 3. Dependence of the real part of the eigenvalues on 
7 and A for a background fluid of matter. Upper and lower 
lines represent respectively solutions — and + in Eq. J2l|). 
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FIG. 4. Dependence of the real part of the eigenvalues on 
7 and A for a background fluid of radiation. Upper and lower 
lines represent respectively solutions — and + in Eq. (pi]). 

For a large class of scalar potentials V(<fi), stability is 



possible for a large number of e-folds allowing a scalar 
field sub-dominance during a long period of time as we 
will see in the applications below. We should point out 
that so far these results are general since we have not 
assumed, yet, any type of potential V(<p). 



V. THE SCALAR FIELD DOMINATED 
SOLUTION 

We now consider the case in which the evolution of the 
scalar field approaches the unit circle x 2 + y 2 = 1, as it 
is the case along region 5 of Fig. 2. Introducing the new 
variables x = XX and y 2 = 1 - X 2 Y 2 , Eqs. (0) give, 



X' = V6A 2 (r - l)X 2 -3X + JUl - X 2 Y 2 ) 2 
-X 2 X[2X 2 + -f(-X 2 + Y 2 )] , 
YY' = V6A 2 (r - \)XY 2 + J-X{1 - X 2 Y 2 ) - 
-(1 - X 2 Y 2 )[2X 2 + 1 {-X 2 + Y 2 )] . 



(23) 



However, the same kind of analysis as for the tracker 
solution becomes here extremely complicated. We take 
then a set of simple reasonable assumptions in order to 
obtain simple and useful results. 

One can see from Fig. |2| that at late times, when the 
scalar field is dominant, y is approximately 1 and x ap- 
proaches zero. This means the scalar potential is over- 
taking the energy density and the potential is very flat. 
In other words, A is getting closer to zero. Let us take 
then A « and A 2 (r — 1) nearly constant. Eq. ( p3| ) then 
reads, 

X 1 = \/6A 2 (r - l)X 2 -3X+ j| , 



(24) 



YY' = V6A 2 (r - l)XY 2 + \j-X 

i -[2X 2 +1 {-X 2 -Y 2 )] 
The system has critical points in, 



*(A)-3| 



y c 2 (A) = i-^. 
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Hence, the scalar field is dominant, fi^ = 1 and 7 ^ 
Aj,/3, where we have defined 



l±^l-4(r-l)A 2 /3 

(r-i)A 



(26) 



for r 7^ 1, and A^ = A otherwise. As before, only the 
minus solution has physical meaning. Expanding A^ we 
can approximate the solution by 



A, 



A 



i + -(r-i)A 2 



(27) 



As we did for the tracker solution, we perturb the so- 
lutions around the critical points to study their stability. 
Expanding Eq. ( j23| ) and using Eq. (E6I) to write T — 1 in 
terms of A and A^, we find the following eigenvalues: 
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m+ = 6 (A* - A)^- + ^{\\ + AA - 67) , 

m_ =3-6 — + -A (3A -2A) 
A,/, I 



(28) 



For r = 1 this expression reduces to the eigenvalues 
found in g], as we would expect. In Fig. || and Fig. ^ we 
show the dependence of the eigenvalues on 7^ and A for 
a background fluid of matter and radiation, respectively. 
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FIG. 5. Dependence of the eigenvalues on 7^ and A for a 
background fluid of matter. Upper and lower lines represent 
respectively solutions m+ and m_ in Eq. (E8|). 
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FIG. 6. Dependence of the eigenvalues on 7^ and A for a 
background fluid of radiation. Upper and lower lines represent 
respectively solutions m+ and m_ in Eq. (E8|). 

For completeness, we indicate that the background 
fluid dominated solution with [x, y, A) = (0,0, A) has 
eigenvalues, 



37 
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771_ = — - 



3(2- 7 ) 



(29) 



and when T = 1, there exist the kinetic dominated solu- 
tions (x, y, A) = (±1, 0, A) with eigenvalues, 



m + = 3 (2 — 7) , m_ 



y|(V6=FA). (30) 



In Table | we give a summary of the properties of these 
critical points. 



VI. APPLICATIONS 



A. Quintessence 



The first estimate for the value of the cosmological 
constant from current particle physics is the Planck scale, 
while the lowest is set by the electroweak scale, namely, 
of order 10 8 GeV . These are tremendously high values 
on cosmological energy scales. So, the simplest solution 
so far has been to assume that by some mechanism, the 
cosmological constant would vanish altogether. 

In the last decade, measurements of the apparent 
magnitude-redshift relation using Snla combined with 
CMB and galaxy clusters and light elements abundances 
measurements, give indications that we are living in an 
accelerating Universe with f2 m attcr ~ 0.3 and J7a ~ 2/3 
0. Then, the discovery that the cosmological constant 
like term is small but non zero, is disturbing. Having 
today a cosmological constant energy density contribu- 
tion, pa = A/87rG, of the same order of magnitude as 
the critical energy density p c ~ 10~ 47 GeV 4 , requires to 
fine tune its initial value to 120 orders of magnitude be- 
low the Planck scale. In order to alleviate this puzzle, 
the idea of "quintessence" was introduced ||. One of 
its versions consistes of an inhomogeneous scalar field Q 
rolling down a potential with attractor solutions. The 
argument is that if the scalar field joins the attractor so- 
lution before the present epoch, information about the 
initial conditions will get lost, which allows a freedom 
in choosing those conditions within 100 orders of mag- 
nitude, thus relieving the fine tuning issue. The initial 
suggestion for such potentials was to use an inverse power 
law form V(Q) ex Q~ a which can be found in models of 
supersymmetric QCD |lj|. Pure exponential potentials 

V tx exp(AQ) also have attractor solutions; however, they 
cannot be used on their own to model quintessence |14| ] 
but interesting modifications have been proposed such 
as V ex exp(— \Q)Q~ a Jig] , sum of pure exponentials 

V tx exp(aQ)+exp(/3Q) ^6j 5 Supergravity inspired mod- 
els with V tx exp(Q 2 )Q- Q @ or V oc exp(Q 4 / 3 )Q 2 / 3 

II- 

In this section we will give the scalar field and its equa- 
tion of state evolution in a background fluid dominated 
setting, for the inverse power law potential. We will then 
turn to the general potential V oc exp(aQ /3 )Q fi . 



Xc 


vl 


Existence 


Stability 


£7^ 


70 




1 

-1 
rw2± 

V2l 






!&(2-7*) 


all 7 and A 
all 7 and A 
all 7 and A 

A 2 > 3 7 ^ 

A^<6 


saddle point for < 7 < 2 

unstable node for A < v6, saddle point for A > y/E 

unstable node for A > — \/Q, saddle point for A < — v6 

stable for Re(m ± ) < (see Eq. @) 

stable for m± < (see Eq. (§|)) 



1 
1 

370 
1 


undefined 

2 

2 
Eq. © 

>i 

3 



TABLE I. The properties of the critical points. A,/, is defined in Eq. (|2 



Consider the inverse power law potential, 
M i+a 



V(Q) 



Q° 



for which the relevant quantities are 



A = 



1 

r- 1 = - 

a 



(31) 



(32) 



We have seen in section Fy] that when the background 
fluid is dominant the equation of state of the scalar field 
can be well approximated by Eq. (pi]). Therefore, in this 
case the equation of state is a constant and is given by 



1Q = 



1 a 



(33) 



Integrating the x equation in Eq. ( jig ) the solution for 
the scalar field reads, 



Q = Q t 



3~t/(2+a) 



(34) 



where Qi can be derived from the y equation assuming 
H 2 « k 2 /3 p l B exp(-37iV), with N = ln(o/Oi), to give 



Qi 



3 kVb 7q(2-7q) 



2 M 4 + a 



-l/(2+a) 



(35) 



We review here some of the properties of power law 
type of potentials ||. For example, Eq. (|33|) still holds 
for negative a provided a < —2. Moreover, we have 
seen that for large A, the solutions are stable provided 
7q < 1 + 7/2, which yields, 



7 + 2 
a > 2 -, a > 0. 



a < 2 



7-2' 

7 + 2 

>-2 : 



a < 0. 



(36) 



For a > the condition is always true, however, for a < 
it imposes the bounds a < —6 and a < — 10 for matter 
and radiation dominated fluids, respectively. For even, 
negative a the field can have early oscillatory behaviour 
with average equation of state given by the virial theorem 



< 1Q >= 1 



a-2' 



(37) 



while energy is continually being converted between ki- 
netic and potential. 

Consider now the more general potential, 



- A/f 4 P Q («Q) 



V(Q) = Af 4 e 



(kQY 



(38) 



we have 



r- 1 



\ = -{al3{nQf + ix){kQ)- 1 
aP{f3 - 1)(kQ) - n 



(aP(KQ)P + M ) 2 



(39) 



and integrating x in Eq. ( |l~6| ) with 7^ given by Eq. (|l£|), 
we find the solution 



a{nQY = $ 



(3 



HkQY 



2\n[al3(KQY + n] - 3~fN , 



(40) 



where the last logarithm term only exists if a ^ 0. $i 
is an integration constant and can be derived from the 
y equation. If we assume kQ 3> [i/af3 then we have a 
nearly constant V — 1 w 0, and we can neglect this con- 
tribution from the equation of state, to yield the simple 
result, 



$,: = In 



2Af 4 



7(2-7) 



(41) 



Since we do not have an explicit expression for (kQ) 13 
we take the ansatz used in [pi, which is a perturbative 
solution, 



cx(kQ)i = 
a{nQt = fy 



$i - 37 TV , 
2-/i 



ln(«Q)Li 



2hx[aP{KQY n _ x +n]-^N, 



(42) 



Substituting this solution back in Eqs. (^9|) and ( p^| ) 
we have now the complete evolution of the scalar field 
equation of state in terms of the background fluid en- 
ergy density and model parameters, only. We can see 
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FIG. 7. Evolution of the equation of state of the scalar 
field in a radiation background. The solid line is the numeri- 
cal evolution and the dashed line our approximate analytical 
solution using Eqs. (ff^), (||) and (jig), 

from Fig. [7| that the second order solution is, in general, 
already a good approximation. 

Note that for a pure exponential ((3 — 1, fi — 0), we 
recover the exact solution 



37 
Q = ft - — N, 



(43) 



where ft — <&i/ na. 

When kQ <C [J./a>f3, a positive definite potential, Eq. 
( pq ) becomes well approximated by a pure power law 
potential, i.e. T — 1 w —1//J,, a ~ and, neglecting the 
second logarithm contribution, we rewrite Eq. (pQ) as, 



kQ 



, a 



3 7 /(2-m) 



(44) 



where ^ = exp($i/(/i— 2)), is estimated, as before, using 
the y equation to yield, 



#< = 



3 Pb_J_JW_ 
2M* fi 2 fi-2 



jfi 



/' 



-l/(2-At) 



(45) 



Despite the above solution being merely an approxima- 
tion, substituting back in Eq. (|39|), the variation of the 
equation of state with time can be well accounted for. 



B. Moduli stabilisation 



effects such as multiple gaugino condensates, which de- 
velop a minimum, and non-perturbative corrections to 
the Kahler potential |lj|. Unfortunately, the scalar po- 
tentials in these theories are exponentially steep, there- 
fore, it is expected that these fields roll past the minimum 
rather than acquiring a vev |20(| . Hence, stabilisation of 
the moduli is one of the most serious questions in string 
theory. Solutions to this problem have been suggested in 
the literature. One possibility is that matter fields other 
than the moduli drive the evolution of the Universe |||| • 
If this is the case, it has been shown that this setting 
opens up a wider region of the parameter space for which 
dynamical stabilisation of the moduli is successful. The 
reason behind this feature is, once again, the existence 
of attractor solutions that slow down the fields by the 
amount needed to trap them when they reach the mini- 
mum. 

Away from the minimum, the evolution of the moduli 
can be written in terms of a canonically normalised field 
a, in a double exponential potential. 

We take the general potential 



V{a)=M 4 eMae pKa ), 



(46) 



where a is an arbitrary real number. Following the same 
line of argument as before, we find the solution for S = 
aexp(/3Kcr) to be (in accordance with B), 



Si 



$, 
$, 



3 7 iV, 
21n(S„_ 1 /a)-3 7 7V, 



and $,; is given by 



$,; = hi 



2 M 4 a 2 j3 2 



(47) 



(48) 



This attractor solution enables us to stabilise the dila- 
ton at the minimum. One sees that the equation of state 
of the scalar field is the same as that of the background 
fluid, apart from a logarithmic deviation and an upper or 
lower shift given by $;. In other words, the ratio between 
kinetic and potential energy of the scalar field is roughly 
a constant and hence prevents the field from rolling past 
the minimum for certain values of the background fluid 
equation of state. 

A last remark is that if a is positive, then A is decreas- 
ing and this becomes a realistic quintessence potential, 
as well. 



In string theory, the modulus and the dilaton (moduli) 
play an important role. They parametrise the structure 
of the compactified manifold and their vacuum expecta- 
tion value (vev) determine the d = 4 value of the gauge 
and gravitational couplings and fix the unification scale 
M GU t- 

Stabilisation of the moduli fields in the d — 4 effective 
theory, has been studied by including nonperturbative 



C. Quintessential inflation 

We have seen in section M that for small enough values 
of A^ the scalar field drives the dynamics of the Universe. 
In this section we present scalar field dominated solutions 
for the same potentials used above. The particular case 
of A^ < 2 is of extreme importance, since it corresponds 
to an inflationary scenario (a tx —p,p — 3p<p > 0). In 



what follows we consider a good approximation, taking 
•V ~ A. 

Conside ring the limit (kQ)^ S> fi/otP, we integrate 
x in Eq.(p5|) to obtain for the potential V(Q) = 
Af 4 exp(a(KQ) /3 )(KQ) A ' the solution, 



( K Qf= [( K Q l ) 2 - l3 + a(3(2-f3)N 
when (3 ^ 2, and 



0/(^-13) 



a(3 



a(3 



-1a(5 



(49) 



(50) 



otherwise (see also |2l]]). The case (kQ) 13 <C [ijaf3 can 
be well approximated by the power law solution, 



{kQY = [(kQ,) 2 - 2fiN] 



0/2 



(51) 



From Eq. (|9|), it should be clear that for (kQ) 13 > 
n/a(3, if (3 < 1, A is increasing and we can imagine a 
scenario in which the Universe is first inflating and then, 
when A is such that A 2 > 2 the universe starts to de- 
celerate. Reheating can happen by gravitational particle 
production by which conventional particles are created 
quantum mechanically from the time varying gravita- 
tional field J22J. Moreover, if we choose negative fi, the 
potential has a minimum with non vanishing vacuum en- 
ergy. Assuming this minimum is the responsible for the 
currently accelerated expansion of the universe, we can 
envisage a scenario of "quintessential inflation" p3[ . 

Taking the number of e-folds from the end of inflation 
to be N — 50 and imposing the spectral index to be 
ns > 0.95, a closer look at this potential reveals that 

/3 < 0.01, (52) 

and for the approximation (kQ) 13 3> —fi/a/3 to be valid, 

-0.1<fj,<0. (53) 

The amplitude of density perturbations ^4s = 2 x 10~ 5 
measured by COBE, then imposes 



a > 293.8, M<10 _x 'GeV 



(54) 



The ratio of tensorial to scalar perturbations is then, 
A^/^lg < 0.17, and the tensorial spectral index nx > 
-0.028. 

An important condition in order not to spoil nucle- 
osynthesis is that the value of the Hubble constant when 
inflation ends must be -ff cn d > 10 8 GeV which is amply 
satisfied in this model as -ff C nd > 10 12 GeV. 

For the double exponential potential, V((j)) = 
M 4 exp(cte° K ^) the field rolls according to 



6 = -In ( e -^ - a(3 2 N) 

k(3 



(55) 



When a is negative the potential offers a plateau for /3(f) < 
0, where inflation can happen. Assuming that particles 



are produced gravitationally as above, to satisfy the value 
of density perturbations we find, M « 10 15 GeV and ng rs 
0.96. Furthermore, for \/3\ > 1, these numbers are nearly 
independent of the precise values of a and (3. We must 
also require, (3 > 0.1 in order not to spoil nucleosynthesis 
predictions. As noted in [J11J, even if the field would reach 
the tracker solution before today, its contribution to the 
total energy density would be ever decreasing. So the 
scalar field can never explain, in this case, the present 
acceleration, unless some other contribution is added to 
the potential to make the scalar field become dominant 
today. 

Although interesting, most of these models of 
quintessential inflation lead to a regime where the ki- 
netic energy is the dominant contribution to the scalar 
field energy density until its value is of order the criti- 
cal energy density today. Therefore, the contribution of 
the scalar field at the time of recombination is negligible 
and the value of the field equation of state today is very 
close to the one of the cosmological constant, 7a = 0. 
Hence, it is not expected that CMB and Snla observa- 
tions will be able to provide direct evidence for these 
scenarios. A second problem associated with these mod- 
els is the overproduction of gravitinos, as the reheating 
temperature Trh is often Trh ~ Hand It 10 9 GeV, lead- 
ing to dangerous cosmological consequences [g4j . And of 
course, a solid theoretical motivation for the suggested 
phenomenological potentials is still to be found. 



VII. SUMMARY 

We have studied the dynamics of a scalar field in a flat 
FRW Universe with a barotropic fluid. We have consid- 
ered the general case of field dependent slope A of the 
scalar potential and analysed the differences, especially 
in the stability conditions, with respect to the constant 
slope case. These results are summarised in table |j. We 
exemplified the tracker and scalar field dominated solu- 
tions by giving explicit calculations of the field evolution 
for different potentials. We have discussed applications of 
these solutions to the issues of quintessence, moduli sta- 
bilisation and quintessential inflation. Finally, we would 
like to emphasise that we have demonstrated here the 
existence and stability of a number of scalar attractor 
solutions which have been previously only assumed in 
the literature. 
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